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Abstract. The designable integr ability (DI) [40] of the variable coefficient nonlinear Schrodinger 
equation (VCNLSE) is first introduced by construction of an explicit transformation which maps 
' VCNLSE to the usual nonlinear Schrodinger equation(NLSE). One novel feature of VCNLSE 

t-H . with DI is that its coefficients can be designed artificially and analytically by using transforma- 

tion. A special example between nonautonomous NLSE and NLSE is given here. Further, the 
I— H ■ optical super-lattice potentials (or periodic potentials) and multi-well potentials are designed, 

which are two kinds of important potential in Bose-Einstein condensation(BEC) and nonlinear 
optical systems. There are two interesting features of the soliton of the VCNLSE indicated by 
the analytic and exact formula. Specifically, its the profile is variable and its trajectory is not 
a straight line when it evolves with time t. 

Keywords: Designable integrability, Variable coefficient Nonlinear Schrodinger equation, 
^ ; Soliton, Optical lattice potential, Double well potential 

in ■ 

CN ; 1. Introduction 

The variable coefficient soliton equations with different formulations have been well studied 
since 1970s, and there are many results on this topic in the literatures. For example, inverse 
scattering method [IH3], symmetry algebra pj], etc. However, the solving and applying of the 
several types of the variable coefficient nonlinear Schrodinger equation(VCNLSE) has been 
^ revived recently in the research of Bose-Einstein condensation(BEC) and nonlinear optics [5HTB] . 
It is well known that the main difficulty in solving process is due to the non-isospectrum 
property in the Lax pair of the integrable VCNLSE. So it is very natural to get the idea that 
one can overcome this difficulty by finding a transformation mapping the VCNLSE to the usual 
nonlinear Schrodinger equation(NLSE) associated with an iso-spectral problem, and thus one 
can get plenty of solutions of the VCNLSE from the abundance of known solutions of the 
NLSE. This idea is pointed out at the added note of the very early work [I] on the NLSE with 
a linear potential of x, which is solved by inverse scattering method. Since then, in order to 
get the explicit solutions of some special cases of VCNLSE with concrete coefficients, there 
are several typical transformations involving dependent variables and independent variables, 
such as lens type transformations [171420] . similarity transformations [714T0 l [T3 |l2TM2"3] and some 
others [2H - [2"9~] . However, for the extensively studied potentials in BEC and nonlinear optical 
systems, i.e. optical-lattice and super-lattice potentials (or periodic potentials) [30H31] and 
multi-well potentials [33 |[35H39] . to our best knowledge, the exact and analytical solutions of 
the VCNLSE such as solitons, have not been reported in the literatures except stationary 
solutions and numerical solutions. In order to show dynamical properties of the solitons in 
one-dimensional BEC and nonlinear optical systems under the control of the two kinds of the 
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external potentials above mentioned, we shall develop a new method to design corresponding 
integrable VCNLSE, and then get its soliton solutions. 

The organization of this paper is as follows. In section 2, a general method which maps 
VCNLSE to NLSE is given with several arbitrary functions. These arbitrary functions provide 
a possibility to design integrable model. In section 3, as a special application of this method, 
the nonautonomous NLSE [8J is designed and mapped explicitly to NLSE. In section 4, two 
NEW integrable models with important physical concerns are designed, and their properties 
are discussed according to analytic solutions given from a single soliton of the NLSE by means 
of our general method. The conclusion will be given in section 5. 



2. General Method 
The investigation object of this article is a very universal VCNLSE in the form of 

d Id 2 

i—ip(x,t) + 2 ^2^( x ' t ) ~ vip(x,t) - g\?p(x,t)\ 2 ?p(x,t) = 0, (1) 

where o = o(x, t), v = v(x, t),g = g(x, t) are three real functions of x and t. This equation is a 
light extension of eq.(l) in reference [13] when 1, which is widely used to characterize one- 
dimensional BEC and nonlinear optical systems under some physical conditions. First focus 
on the integrability of eq.([T]) by means of looking for a direct relationship between VCNLSE 
eq.([T]) and usual nonlinear Schrodinger equation (NLSE) eq.([2]). Then, find a transformation 
mapping eq.(0Q) to the usual NLSE, 

■ dq I d 2 q 2 

t df + 2dX 2+elqlq = ^ (2) 

where e = ±1, q = q(X,T). Meanwhile, coefficients o,v,g are given analytically by this 
transformation. Therefore, one advantage of this transformation is to solve VCNLSE by using 
all known solutions of NLSE as we discussed in the above section. 
To this purpose, a trial transformation 

tl>(x,t) = q(X,I)p(x,t)J«**) (3) 

is introduced with X = X(x,t),T = T(t). So the central task is to determine the concrete ex- 
pressions of real smooth functions {o, v, g, X,p,<f>} by requesting q(X, T) to satisfy the standard 
NLSE eq.(T5]). By substituting the transformation eq.([3]) into eq.([T]), and setting 

o(x, t) = — g(x, t) = (4) 

without loss of generality of the transformation, then it becomes 

• , 1 . | i2 , 1 (**3 + h)q . 1 (iki + k 2 )q x n ,^ 

"> t + 2 qxx +eMq+ 2im^ + 2 w = °- (5) 

dq dT 

Note qx denotes — — , T t denotes — — and so on. Here fcj, i — 1, 2, 3, 4, are given by 
oX at 

h = 2pX t X 2 + 2pT t <p x X x , k 2 = T t (2p x X x + pX xx ) , 

k 3 = 2p t X 2 + T t (2p x (/) x + p4> xx ) , k 4 = T t p xx - 2pcj) t X 2 - T t p4> 2 x - 2vpX 2 . 
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Obviously, let k\ = k 2 = k 3 = fc 4 = 0, then eq.flS]) becomes the standard NLSE. Moreover, we 
get 

X = X(x,t) = J F{x)f 1 (t)dx + f 3 {t), T — T(t), (6) 

P = p(x,t) = ^^= (7) 
y/F{x)jx{t) 

* M *\ f {jF(x)f lt dx + f 3t )F(x)f 1 (t) ^ 

<f> = <f>{x,t) = - J — dx + f 2 (t), (8) 

. l-3T t F 2 + 2T t F xx F + 8f 1 (t) 2 <f )t F* + 4T t (<p x ) 2 F 2 

v{x > t] = "s JWfW ' (9) 

from k\ = k 2 = k 3 = k± = by a tedious calculation. Furthermore, taking X, p, T back into 
eq.flU), then 

/x T t , eF(x)T t , , 

o(x,t)= \ g(x,t) = — io 

So the transformation in eq.()3]) indeed maps the VCNLSE to the NLSE as we wanted, which is 
determined by eq.([6]) to eq. ffTOj) through five real arbitrary functions fi(t), f2(t), f3(t),T(t), F(x), 
F(x)f\(t) > 0. At last, we would like to point out that we may set X = X(x, t) and T = T(x, t) 

d 2 

in transformation eq.© for better universality. However, to eliminate the term ( q(X, T)) 

d d 

in the transformed equation, we have to ask (—X) = or (— — T) = 0. So we choose directly 

ox ox 

T = T(t) for simplicity in eq.©- The other choice X = X(t),T = T(x,t) will be given in a 
separate paper. 

Thus we can design o, the external potential v and interaction nonlinearity g according to 
different physical considerations by means of the selections of the arbitrary functions above 
mentioned, such that the integrability is guaranteed. Therefore, we call that the VCNLSE 
possesses the designable integrability(DI) [40J, which originates from the rigid integrability [ID] 
of the NLSE and the transformation eq.Q. 

3. NONAUTONOMOUS NONLINEAR SCHRODINGER EQUATION 

By comparing with the known results on the connection between VCNLSE and NLSE men- 
tioned at the first paragraph, our result is more universal because the research equation eq.flTJ) 
and transformation eq.OH]) is more general. Moreover, the integrable conditions [7]l8"| l2"oT - [2"8"] of 
the coefficients in the VCNLSE disappear in our method. Of course, to guarantee the integra- 
bility of VCNLSE, these coefficients can not be arbitrary functions. Actually, these conditions 
are satisfied automatically by analytical expressions of coefficients o, v and g. 

To show this point clearly, we would like to demonstrate that the nonautonomous NLS j8] has 
designable integrability. In other words, we shall show how to choose functions fi(t), f2{t), f3(t), 
T(t),F(x) such that we can get nonautonomous NLSE [8] 

z^ + ^^ + eR(t)\Q\ 2 Q-2a(t)xQ-^x 2 Q = 0, e = ±1, (11) 
with a integrability condition 

- = Tt ln c « + R ^W) - S lnD K laR ^ < 12 > 



from VCNLSE, and nonautonomous NLSE can be transformed to the usual NLSE by pre- 
vious transformation eq.()H]). So by comparing with nonautonomous NLSE, then g(x,t) = 



—eR(t),o(x,t) = D(t),v(x,t) = v 2 x 2 + Vix + v Q ,v 2 = — - — , t>i = 2a(t),v = 0. To get this 
integrable model, setting F(x) = F , and using o(x,t) and g(x,t) we get 

= ft- (13) 



T(t)= I -xn^dt + To- (14) 



R{tf 

D{t)F$ 

Here Fq is a real constant, To is a integral constant. Taking F(x) = Fq into the expression of 
v(x,t) in eq.(jH]), it infers that v(x,t) is a second order polynomial of x. Then 

f2(t) = -i / 1 } ^ cK + c 20 (15) 
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and 

D(t)a(t) 
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/s(«) = y - — R{ % D{ t) — : — dt + Cso2 (16) 

are given from eq.()9]) with the help of V\ = 2a(t),v = 0. Here c 2 o, c 30 i, c 30 2 are integral 
constants. Moreover, taking v 2 = ^{t) 2 in v(x,t) infers 

1 R t D t R 2 1 D tt 1 R tt 1 D 2 1 2 

2 i?D 2 i? 2 D 2 D 2 2 RD 2 2 ' K J 

which is equivalent to eq.( !T2"]) . Further, taking F(x) = F and fi(t) given by eq. (fT3l) into eq.([7j), 
then 



'-M-^WY (18) 

According to eq.© and eq.flHJ), and using known functions fi(t), f2(t), f 3 (t),T(t), F(x), then 
X= Fof^x + Mt), (19) 
0= 6 2 (t)a; 2 + 6i(t)a; + 6o(t), (20) 
u _ l^-ADt) (21) 



2 



2 



6 '= > dW ~ (22 » 

6o = AW = "J _/ W D(t) dt + C2 „. (23) 



We have verified that the transformation 



Q[x, t) = q{X, T)p(t)exp(i<P) (24) 
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given by eq. f lT4"j) ,eq. (ITS]) .eq. (IT9"]) and eq.t lZO]) indeed maps nonautonomous NLSE eq. ffTTj) to the 
usual NLSE. It is trivial to find that transformation eq. (}24|) is equivalent to the result in 
Ref. [26]. 

As the end of this section, we would like to show several examples of nonautonomous NLSE 
eq. tfTT]) . and integrable condition Eq.( fl2"]) held for them. Their solutions can be obtained from 
known solutions of the NLSE by transformation Eq. (I24j) . In particular, all of the following 
equations has DI property. 

Example 1 [15]: Let D(t) = 2,R(t) = 2g e xt , f2(i) 2 = — ^-,a(t) = 0, then the nonautonomous 
NLSE reduces to Eq.(l) of reference [15] with a(t) = goe xt . 

Example 2 [41J : Let D(t) = l,R{t) = l,Q{t) 2 = Q,a(t) = \{B + Csin(itrf)), then the 
nonautonomous NLSE reduces to Eq.(2) of reference |41j . 

Example 3 02]: Let D(t) = 2,R(t) = 1 + tanh(wt), Q(t) 2 = (tanh(u;t) - l)w 2 ,a(t) = 0, then 
the nonautonomous NLSE reduces to Eq.(l) of reference [12] with g(t) = 1 + tanh(iot) and 
k(t) = 1 -tanh(wt). 

Example 4 HUDS]: Let D(t) = l,R(t) = 1 + m sin(wt), Q{t) 2 = 2k(t),a(t) = 0,then 
the nonautonomous NLSE reduces to Eq.(2) of reference [19] with g(t) = 1 + msin(wt) and 

1 mw 2 (2m cos(wt) 2 + sin(wt) + m sin(wt) 2 ) 

2 1 + 2msm(wt) + m 2 sm.{wt) 2 

Example 5 [43]: Let D(t) = 1, R{t) = 1, Vt{t) 2 = 0, a{t) = -f{t)/2 ,then the nonautonomous 
NLSE reduces to Eq.(3) of reference [13] with rj = — 1 and f(t) = b\ + / cos ut given by eq.(30) 
of reference [43J. 



4. New Examples 

To further illustrate the wide applicability of our methodology, and motivated by the ex- 
treme importance of the external potentials in the BEC and nonlinear optics systems, two 
NEW integrable models: integrable VCNLSE with optical super-lattice potentials (or periodic 
potentials) and multi-well potentials are designed respectively. In the following examples, we 
shall set = l,/ 2 (t) = l,/ 3 (t) = l,T(t) = t,e = l. 

New Example 1: Optical super-lattice potentials According to transformation eq.([3]), we 

3 1 3 

have three potentials:l)F(x) = l/(cos(:r)+3/2), v(x, t) = cos(2x) — cos(a;); 2)F(x) = 

3 93 11 39 

1/ (cos(2a;) + cos(x) + 21/10), v(x, t) = —2 cos(x) 4 cos(x) 3 cos(x) 2 — — cos(x) + — ; 

2 40 40 40 

4 10 

3)F(x) = 3/(cos(3x) + cos(2x) + 5 cos(x) + 10), v (x, t) = —2 cos(x) 6 — - cos(a;) 5 — — cos(a;) 4 — 
80 25 11 17 

— cos(x) 3 — — cos(a;) 2 + — cos(x) + — from eq.([n]). They are real periodic functions of x with 

period 2tt, and have one peak, two or three peaks over intervals of length 2n respectively. The 
profiles of them are plotted in Figure 1 from left to right in order. 

New Example 2: Multi-well potentials According to transformation eq.OH]), we have got 

three types of multi-well potentials: 1) Let a > 0, then Fix) = ; —, — ; — rr 

yi 1 ! y ' 200cos(arccot((x/20) 2 + a)) 

gives symmetric double well potentials, which are plotted in Figure 2 (left) for a = 0.001,0.2 

0C CC tjC 

and 1; 2)Let a > 0, then F(x) = ((ttt:) 4 + a ("P7:) 2 + (37:) + 4)/400 gives non-symmetric double 

oU OU oU 

well potentials, which are plotted in Figure 2(middle) for a = 1/15,2 and 7. Note that the 
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Figure 1. The designed optical super-lattice potentials(New Example 1) with 
period 2tt. 
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last one is a single well potential. 3) Fix) = ; , ...... ; — - , ,„ . . 

& 1 ' K J 15 2cos(4sech(§)) + 1.3cos(sech(§)) + 10 

triple well potential, which is plotted in Figure 2(right). 
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Figure 2. The designed multi-well potentials(New Example 2). Left 
panel(the symmetric double well potentials): From the thinnest (magenta) to 
the thickest (green) curves it denotes potential generated by F(x) with a = 
0.001,0.2,1, respectively. Middle panel(the non-symmetric double well poten- 
tials): From the thinnest (magenta) to the thickest (green) curves it denotes po- 
tential generated by F(x) with a = 1/15, 2, 7,respectively. Right panel:the triple 
well potential. 



Note that the associated other two coefficients o(x,t) and g(x,t) are given simultaneously by 
means of transformation eq.flSJ, such that VCNLSEs associated with these designed coefficients 
o, g and v are integrable systems. To save space, we do not write out them here. Moreover, 
double well potential can be achieved by taking F(x) = 0.35x 2 + 5sech 2 (0.27x). Additionally, 
by setting arbitrary functions in transformation eq.([3]), the interested readers can design dif- 
ferent integrable VCNLSE to establish mathematical model equation of physical systems. For 
instance, time-dependent optical super-lattice potential jH] can also be designed by setting 
T = f f(t)dt and F(x) = l/(acos(x) + b) in v(x, t). Here b > a > 0. 

Furthermore, the transformation eq.Q provides an efficient way to construct exact and 
analytic solutions if)(x, t) of the designable VCNLSE from known solutions q(X, T) of the NLSE, 
such that we can explore the dynamical evolution of solitons of the VCNLSE conveniently. For 
example, setting a usual single soliton q(X, T) as 



\q(x,T)\ 



2\/2~M 



cosh(2r7iV2X + 877i£ x T) 



(25) 



which is given by eq.(6) of reference [24]. Obviously, the profile of the usual soliton q(X,T) 
on the plane of (X,T) is invariant when soliton evolves with time t. A single soliton of the 
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VCNLSE, 

\^(x,t)\ = \p(x,t)UX,T)\ = \p(x,t)\ 2 ^ (26) 

cosn(2?7iv2A + 8^i4iT) 

is obtained by using transformation eq.((3]). Note that the profile of the soliton \i[)(x,t)\ of 
the VCNLSE is not preserved when it evolves with time t because the amplitude \p(x, t)\ is a 
function of x and t and the trajectory is a curve x = x(t) on the plane of (x,t), which is defined 
implicitly by 

V2X(x,t)+4^T(t) = 0. (27) 

This shows that the profile of the soliton of the VCNLSE is designable by using different p(x, t), 
X(x, t) and T(t), which can be realized by choosing different arbitrary functions fi(t), fs(t), F(x) 
and T(t) in transformation eq.([3]). In particular, as we shall show in the following example, 
the amplitude of ip(x,t) is dependent of t even if the p = p(x) is x-dependent only, because 
p = p(x(t)) when soliton moves along the trajectory x = x{t). 

To further illustrate the property of dynamical evolution of the soliton \ip(x,t)\ of the VC- 
NLSE, one example is given here. According to the above formula eq. (l2"6"|) . the solution of the 
case 1) in New Example 1 is a deformed single soliton, i.e., 

I#M)|= , , 2MI ■ (28) 



^=rZIcosh(2V2, 1 X + 8 fl ^ 



~ 4 1 x 

which is plotted in Figure 3 with £i = 0.1 and r/i = —0.1. Here X = -^(arctan(— F tan(— )) + 
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FIGURE 3. The dynamical evolution on (x — t) plane of a single soliton of the VCNLSE with 
an optical super-lattice potential designed by case 1) of New Examplel. 

X _|_ 7T 

[- — 1 7TJ + 1, and \x] is the greatest integer less than or equal to x. Particularly, X is a 

7T L J 

monotonically increasing and continuous function of x, which is obtained from eq.([6]) by adding 
the floor function [x] . Figure 3 shows intuitively two interesting features of the soliton of the 
VCNLSE: 1) the profile is variable and 2) the trajectory is not a straight line when it evolves 
with time t, as we pointed in the above paragraph. This is also supported visibly by Fig.4 for 
profiles of different time t and Fig 5 for the trajectory on the plane of (x,t), which is given by 
V2X + OAt = 0. 

Note that there are many other kinds of solutions of NLSE, such as dark soliton, periodic 
solution, position, negation, complexiton, etc., which can be used to generate the solutions of the 
VCNLSE. In ourexamples, merely the single bright soliton of the NLSE is used. Especially, from 
the "seed" of the multi-soliton solutions of NLSE, the multi-soliton solutions of the VCNLSE 
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FIGURE 4. The one-dimensional profiles of single soliton in Fig. 3 at different times. From 
the right (magenta) to the left(green) curves it denotes soliton at t = —55,0,55, respectively. 
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FIGURE 5. Trajectory of the single soliton in Fig 3. 




can be obtained, and then its interaction properties are discussed. However, it is very difficult to 
find multi-soliton solutions of the VCNLSE from the view of non-isospectral integrable system. 

5. Conclusions 

In conclusion, a new concept "designable integrability " jlQ] of the VCNLSE has been devel- 
oped by a novel transformation eq.Q. The novel characteristic of VCNLSE with DI is that the 
variable coefficients with physical meaning including o(x, t), v(x, t) and g(x,t) can be designed 
artificially and analytically according to different physical considerations. Correspondingly, 
the profile of its solutions including soliton and others are also tunable intentionally by us- 
ing different p(x,t), X(x,t) and T(t), which can be realized by choosing different arbitrary 
functions. The nonautonomous NLSE [8J is re-obtained as a special design through choosing 
F(x) = Fq (constant) . Furthermore, two kinds of NEW VCNLSE with optical super-lattice 
potentials and multi-well potentials have been designed respectively, and two interesting fea- 
tures of the soliton of the VCNLSE are shown by formulas and figures. The results in this 
paper show that DI and some unusual behaviors of solitons for the VCNLSE originate from the 
usual NLSE and the transformation eq.(j3]). The methodology of studying variable coefficients 
partial differential equations with DI can be extended to many other cases including other 
(l+l)-dimensional, even higher dimensional integrable systems, multi-component systems and 
discrete systems. In particular, the method to design the solvable o, v and g of the VCNLSE 
is expected optimistically to be used by theoretical and experimental researchers. 

As the end of this paper, the method applied to a more general VCNLSE 

t) + l ^){x, t) - vi/j(x, t) - g\iP(x, t) | 2 ^(x, t) - iTiP(x, t) = 0, (29) 
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which is a minor extension of modeling equations of BECs and nonlinear optics in recent many 
works, is shown here. In order to guarantee integrability of eq. (|29p . Luo et al have shown that 
o, g, T must be time-independent, and v (x, t) must be a second order polynomial of x and satisfy 
a integrable condition (eq.(22) of reference [28]). However, their integrable conditions are too 
strong because of the restriction of the Painleve analysis, and a very general forms of o,g,v , 
which are real functions of x and t, has been achieved by a similar transformation of eq.(j3]). 
This provide us more possibility for soliton control, which will be given in a separate paper. 
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